/// - ? o 

// O / 


An Investigation of the Largest Flares 
in Active Cool Star Binaries with ALEXIS 

NASA P.O. No. S- 68399-F 
Final Report 
Submitted To: 

National Aeronautics and Space Administration 
Goddard Space Flight Center 
Greenbelt MD 20771 


Principal Investigator: 

Robert A. Stern 
Dept Hl-12 Bldg 252 
Solar and Astrophysics Laboratory 
Lockheed Martin ATC 
3251 Hanover St. 

Palo Alto, CA 94304 
(415) 424-3272 


June 3, 1998 


1. Summary 

We proposed to search for high-temperature, flare-produced Fe XXIII line emission from 
active cool star binary systems using the ALEXIS all-sky survey. Previous X-ray tran- 
sient searches with ARIEL V and HEAO-1, and subsequent shorter duration monitoring 
with the GINGA and EXOSAT satellites demonstrated that active binaries can produce 
large (EM ~ io 56-56 cm -3 ) X-ray flares lasting several hours or longer. Hot plasma from 
these flares at temperatures of 10 7 K or more should produce Fe XXIII line emission at 
A = 132.8A, very near the peak response of ALEXIS telescopes 1A and 2A. Our primary 
goals were to estimate flare frequency for the largest flares in the active binary systems, 
and, if the data permitted, to derive a distribution of flare energy vs. frequency for the 
sample as a whole. 

After a long delay due to the initial problems with the ALEXIS attitude control, 
the heroic efforts on the part of the ALEXIS satellite team enabled us to carry out this 
survey. However, the combination of the higher than expected and variable background 
in the ALEXIS detectors, and the lower throughput of the ALEXIS telescopes resulted 
in no convincing detections of large flares from the active binary systems. In addition, 
vignetting-corrected effective exposure times from the ALEXIS aspect solution were not 
available prior to the end of this contract; therefore, we were unable to convert upper 
limits measured in ALEXIS counts to the equivalent L buy- 


2. Technical Progress 

2.1 Background 

X-ray emission from RS CVn and related active binary systems is commonly known to 
be much stronger than from typical single star coronae (Pallavicini 1989; also see the 
Strassmeier et al. 1993catalog and references therein). These same systems also produce 
some of the largest known stellar flares in terms of X-ray luminosity or emission measure. 
Stem (1991) reviewed recent studies of X-ray flaring in such systems: Figure 1 from this 
paper is reproduced below. 

A quick inspection of Figure 1 reveals several important facts: (1) even the largest 
solar flares pale in comparison to the largest RS CVn flares, (2) the classical “flare stars” 
also fail, by orders of magnitude, to reproduce the emission measures of the RS CVn 
flares, and (3) at least in terms of X-ray continuum temperature, the largest flares also 
seem to have the highest temperatures. 

There were a limited number of RS CVn flares seen with Einstein, including one 
on <t 3 CrB (Agrawal et al. 1981), and another on the RS CVn-like Hyades binary HD 
27130 (Stem et al. 1983). However, the most detailed studies of stellar flare properties 
in active systems were made with EXOSAT, and, more recently, GINGA. Flares in the 
active binary systems <r 2 CrB , II Peg, TY Pyx, and Algol were seen by EXOSAT (see 
review by Pallavicini 1989). Some of the brightest X-ray flares were observed with the 
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Figure 1: (a) Emission measure {n\V ) vs L x for a variety of solar amd stellar flakes 
(G=GINGA data). Approximate ramges for dMe flares aind a group of lairge RS CVn 
flaires seen by Ariel V are shown. The L x vs EM relation estimated by Tanaka(1987) 
for solar flaires is shown ais a dotted line, (b) Si mil ar plot for flare temperature. The 
approximate relations for “high T” and “low T” solar flatres suggested by Tanaka (1987) 
are shown as dotted lines. In both plots, shaded region indicates approximate ramge of 
solair data. 

GINGA satellite in UX Ari (Tsuru et al. 1989), II Peg (Doyle et al. 1991) and Algol 
(Stern et al. , 1992). A summairy of large flaires seen in RS CVn systems is given in Table 
1 . 

Despite their limi ted number, the EXOSAT amd GINGA observations yield impor- 
taint constraints on flare physicail parameters, including temperatures, cooling times, 
densities, amd size scades. There is, however, considerably less knowledge regarding flaire 
frequency in such systems. This is simply because the vast majority of information we 
have regarding stellair X-ray flaures was obtained with pointed instruments on Einstein, 
EXOSAT, and Ginga. Yet events with am emission measure (EM=n*7 ) > 10 64 cm -3 
camnot be too rare; otherwise, the success rate for detecting flares in the very limited 
number of several-day pointed observations would be nil. As in the case of solair flaires, 
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Table 1: Largest Flares Reported for RS CVn Systems and Algol 


System 

L*° 

Reference 

<r 2 CrB (EXOSAT) 

« 5 

van den Oord et al. (1988) 

HR 1099 (Copernicus) 

22 

White et al. (1978) 

HR 1099 (Ariel V) 

600 

Pye and McHardy (1983) 

UX Ari (Ariel V) 

600 

Pye and McHardy (1983) 

UX Ari (em GINGA) 

120 

Tsuru et al. (1989) 

AR Psc (HEAO-1) 

250 

Ambruster and Wood (1986) 

II Peg (Ariel V) 

150 

Pye and McHardy (1983) 

II Peg (em GINGA) 

30 

Doyle et al. (1991) 

o Gem (Ariel V) 

400 

Pye and McHardy (1983) 

Algol (EXOSAT) 

« 10 

White et al. (1986) 

Algol (em GINGA) 

10 

Stern et al. (1990) 


“lO 30 ergs" 1 , 2-10 keV 

there is likely to be a “log N-log S” type relation for flare frequency vs energy: on the 
Sim, the number of X-ray bursts above a given threshold scales roughly as S _1 , where 
S is the peak energy flux (e.g., Drake 1971). At this time, there are simply too few 
observations to determine such a relationship for active binary X-ray flares; what little 
data exist are from X-ray transient surveys. 

Although there have been many studies of X-ray transients, the two most sensitive, 
and therefore most useful in estimating flare rates, were the Ariel V SSI Jill sky transient 
survey (Pye and McHardy 1983) and the HEAO-1 A-l all sky transient survey (Am- 
bruster and Wood 1986). The HEAO-1 survey had a sensitivity limi t of ss 10" 10 erg 
cm -2 s -1 in the 0.5-20 keV range, while the Ariel V sensitivity limi t (for a single orbit) 
was ~ 3-4 x lO -10 erg cm -2 s' 1 in the 2-18 keV band. At a distance of 50 pc, these 
sensitivity limits correspond to 3 x 10 31 to > 10 32 erg s -1 . Both studies were carried 
out with scanning collimated proportional counters: the HEAO-1 instrument had a 1° 
x 4° FOV, and the Ariel V SSI a 0.75° x 10.6° FOV. Both missions operated in a sky 
scanning mode, with the long axis of the instrument field of view either perpendicular 
(HEAO-1) or highly inclined (Ariel V: 65°) to the spin plane. 

Both experiments have detected the largest stellar flares observed to date (Table 1). 
A typical position in the sky was in the detector field of view only a limited number 
of days per year (up to 20 in the case of Ariel V, only 8 or so for the 7 months of the 
HEAO-1 A-l survey). In spite of this, at least one of the HEAO-1 and six of the Ariel 
V transients were identified with RS CVn flares with L a = 10 32-33 erg s -1 (Ambruster 
and Wood 1986, Pye and McHardy 1983). In the temperature range 20-60 x 10 fl K, 
typical for the X-ray emitting plasma in such flares, the 2-10 keV band emissivity is ~ 5 
x 10~ 24 to 10~ 23 erg cm 3 s -1 (Mewe, Gronenschild, and van den Oord 1985). Hence the 


3 


EM for these flares is « lo 65-56 cm -3 (see Figure 1). Pye and McHardy (1983) also note 
that the typical duration of the Ariel V flares is ~ 3 orbits, or roughly 5 hours, with the 
notable exception of a 12 hour flare on cr Gem. A more recent example of a long flare is 
the moderately large (EM~10 54 cm -3 ) Algol event seen by GINGA, with an exponential 
decay time of ~ 10 hours (Stern et al. 1990), as shown in Figure 2. 
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Figure 2: Temperature (solid line) and emission measure (dotted line) vs. time during 
Algol flare seen by GINGA (From Stern et al. 1990). 

All the above results lead us to the conclusion that large flares are occurring with 
enough regularity on active binaries to classify them as true “flare stars,” but on a much 
grander scale than the classical dMe systems. The long flare duration and huge emission 
measures also provide evidence that detection of similar or even fainter flares will be 
possible using the ALEXIS multilayer telescopes. Plasma at temperatures in the 10 7 
K and higher range produce line emission, not only at X-ray wavelengths, but also in 
the EUV. In particular, the Fe XXIII line at 132.8 A has a peak emissivity near 10 7,2 
K (Mewe, Gronenschild, and van den Oord 1985, Landini and Monsignori-Fossi 1990). 
Hence the cooling flare plasma from large X-ray flares should be seen as an increase in 
count rate in the 93 eV (130A) band of ALEXIS. 


4 




2.2 ALEXIS Launch Problems and Recovery 

Although this program was first proposed in 1992, and ALEXIS was launched in 1993, 
funding was delayed until March 1996 because of ALEXIS launch problems and the usual 
contractual delays. At some point during the ALEXIS launch, one of the brackets holding 
a solax panel (and magnetometer) to the ALEXIS spacecraft was broken. This prevented 
the satellite from properly acquiring the Sun during its start-up phase, and almost (but 
not quite) draining the batteries. As the ALEXIS orbit processed, the position of the 
solar arrays relative to the Sun gradually changed, allowing initial acquisition of the 
satellite roughly two months sifter launch. 

Remarkably, the solar panel with the broken attachment is still attached to the 
spacecraft via its wire bundle; however, this loose attachment has resulted in a significant 
“wobble” to the spin stabilized satellite and the need for constant ground adjustment 
of the spin axis orientation(Roussel-Dupre et. et. 1997). Hence significant effort was 
required to obtain attitude solutions. These efforts were successful, and ALEXIS is able 
to obtain positions of several bright steady EUV sources (such as HZ43) to ~ 0.2 deg. 

2.3 In-Orbit Performance 

The in-orbit performance is not as good as expected pre-launch because of a number of 
factors: 

(1) Telescopes 1A and 2A have roughly 60% of the response expected from pre-flight 
models (however, much of this is due to a narrower bandpass, which is only of 
significance for off-peak wavelengths) 

(2) ALEXIS discovered a significant in-orbit background correlated with the ram angle 
of the satellite in orbit which results in background levels thousands of times greater 
than predicted. However, this background is present only during roughly 50% of 
each orbit; thus “good times” can be selected from the data, albeit with the effective 
exposure times per orbit reduced significant y, thus decreasing the sensitivity. 

The sensitivity of the individual ALEXIS telescopes to coronal plasma is for am emis- 
sion meaisure of 10 52 cm -3 at 10 pc is shown in Figure 3. 


5 



Figure 3: ALEXIS Sensitivity 


2.4 Data Processing Approach 

The data processing approach for this project differed somewhat from that adapted for 
the “standard” ALEXIS transient searches. Specifically, (1) the target list was known 
a priori , and (2) the time scales expected for the peak emission of large stellar flares 
were expected to be ~ one to several ALEXIS orbits. Thus summing observations over 
a day or more (used in the case of some transient searches) would actually reduce the 
signal-to-noise of the observation. During a meeting at Los Alamos in April 1996, we 


6 






(J. Bloch and the LANL/ALEXIS team, X. Wu of the UCB/ALEXIS team and myself) 
agreed upon a modified strategy to detect stellar flares: 

• For each target, a “pigeonhole” of aspect-corrected photon events was generated 
in sky coordinates with a 5° radius and 0.25° pixels. 

• Each “pigeonhole” consisted of data taken in each of the 6 telescopes from Decem- 
ber 1994 to April 1996. 

• The event list was corrected to include only times of low anomalous background. 

• For each “pigeonhole”, a source region of 30 arc-min radius and a background 
annulus from 50-200 arc-min radii were chosen to test for the presence of a source 
using the method of Lampton (1994). 

• The Lampton “score” for a single orbits worth of data and sums up to the 8 prior 
orbits were calculated and tabulated. 

• As a statistical test of the technique, artificial data sets were generated using the 
same photon data but the aspect information rotated by 180° in the sky. 

• As a test of the “efficiency” of the search techniques, we included in the target list a 
few sources known to have been previously detected in quiescence by the standard 
ALEXIS source detection algorithm: e.g. HZ43. 

2.5 Target List 

The list of targets used to create pigeonholes is given in Appendix A. The list consist of 
essentially all the objects in the Strassmeier et al. ( 1993) catalog. In addition, we included 
in the list the EUVE bright sources mentioned above, four nearby Algol systems, and 
(for another project headed by S. Cully using the same data set), a list of nearby flare 
stars and two cataclysmic variables. Results from the latter sources were ignored in the 
following analysis. 

2.6 Results 

Comparison of Original and “Rotated” Data Sets 

The bright EUV sources put in our target list as test cases were successfully detected in 
quiescence by our transient search algorithm. At first, we also thought that one or more 
transients had been detected at a high level of significance: however, some detective work 
by Jeff Bloch demonstrated that these obvious detections were, in fact, due to reflected 
solar EUV radiation from the Moon. 

After removal of the “test” detections from the EUV bright source list and the 
serendipitous lunar EUV detections, we examined the histogram of Lampton scores for 
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a portion of the data set. This consisted of roughly 120,000 single orbits of data, divided 
unequally among the 6 telescopes. We added together the telescope combinations with 
the same bandpasses (1A+2A, 1B+3A, and 2B+3B), and compared the actual data 
set with one whose aspect information had been artificially rotated by 180° about the 
telescope spin axis, thus, in principal providing a control sample to check the detection 
statistics. For each single orbit, we also added up the previous 2-9 orbits with acceptable 
data for each “pigeonhole”, in an attempt to check our ability to detect longer transients. 
The results of this series of tests is presented graphically in Appendix A, Figures 1-9. 
Note that, in essentially all cases, the original and 180° rotated data sets produce the 
same distribution, with no Lampton Scores above 6.0. Recall that we have removed all 
“bright” steady sources such as HZ43 and the few objects known to be contaminated by 
the Moon. Thus, there is no strong evidence for a non-statistical distribution of transient 
sources with this portion of the data set. 

Histograms of All Processed Data 

Continuing in the above vein, we processed all the data then available. This represented 
approximately 16 months worth of data, or ~ 1.6xl0 8 telescope-orbit combinations. 
Hence, any detections significantly above a Lampton Score of about 6.0 for single orbit 
data should be clear transients. Once again, we summed up to 9 previous orbits to 
search for longer time scale events: note that the latter searches are not independent 
of the single orbit seraches or of each other. The results of these searches is shown 
in Figures 1-9 of Appendix B. In this case, we detected a number of events above the 
Lampton Score = 6.0 value; for single-orbit seraches, we would expect about 1, so, at 
first glance, the results look promising. However, a number of these “detections” me 
also lunar detections for objects such as V471 Tau and BD 17° +703. After removal 
of these, we detect about 5 objects above the 6.0 threshold. Since this includes all the 
multiple orbit searches, such a number is not particularly significant. Thus we are forced 
to conclude that ALEXIS has not detected any coronal transients with a high level of 
confidence. 

Conversion to Flux Upper Limits 

Given the above null result, we still could convert our data to provide upper flux limits 
in each of the telescope bands for our search. However, the software and processing 
required to determine effective exposure time and vignetting corrections to convert the 
counts into erg cm -2 s _1 was still in the process of being developed at the UCB ALEXIS 
Data Center at the time of completion of this contract. The data are still available, so 
that this task could be completed at a later date. 
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2.7 Summary and Conclusions 

We searched « 16 months of ALEXIS data for coronal flare-like transients, but obtained 
no unambiguous detections. The quantitative conversion of our results into upper limits 
awaits the completion of ALEXIS effective exposure software. 
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Appendix A 


Comparison Histograms: Original 
and 180 ° Data Sets 


Number of Hits Number of Hits Number of Hits 





Figure 1: Comparison of Original and 180° Data: 1 Orbits Summed 
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Figure 2: Comparison of Original and 180° Data: 2 Orbits Summed 
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Figure 3: Comparison of Original and 180° Data: 3 Orbits Summed 
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Figure 4: Comparison of Original and 180° Data: 4 Orbits Summed 
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Figure 6: Comparison of Original and 180° Data: 6 Orbits Summed 
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Figure 7: Comparison of Original and 180° Data: 7 Orbits Summed 
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Figure 8: Comparison of Original and 180° Data: 8 Orbits S umm ed 
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Appendix B 

Histograms of Summed Data Set 
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Figure 6: Histogram of Lampton Scores for All Data: 6 Orbits Summed 
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Appendix C 
Target List 


# RA (2000) DEC NAME 


107.515 

38 . 543 

G1 268 

169.545 

31.529 

xi UMa B 

339 . 678 

-20.623 

FK Aqr 

283.863 

8.404 

VI 285 Aql 

1.422 

45.810 

ADS 48A 

223 . 349 

19.153 

HR 5553 

96 . 543 

18.757 

0U Gem 

38 . 591 

-44.798 

CC Eri 

79.173 

45.998 

alpha Aur 

113.658 

31.870 

YY Gem 

315.012 

40 . 074 

V1396 Cyg 

298 . 573 

-23.941 

HR 7578 

248.601 

57 . 149 

CM Dra 

278.483 

51.719 

BY Dra 

69.198 

27 . 134 

V833 Tau 

20.737 

7.419 

AR Psc 

248.680 

-4.219 

G1 629. 2A 

8.812 

-3.593 

13 Cet A 

243.669 

33.858 

sigma 2 CrB 

42.182 

31.115 

VY Ari 

347 . 489 

47.958 

KZ And 

174.320 

47.453 

DF UMa 

354.391 

46.458 

lambda And 

283 . 972 

23 . 557 

V775 Her 

10.697 

35 . 547 

FF And 


D(PC) 

PER (D) 

SP TYPE 

5 

10.428 

dM5e/dM5e 

8 

3.980 

G5V 

8 

4.083 

dM2e/dM3e 

11 

10.319 

M3.5Ve/M3.5Ve 

11 

0.000 

dK6 

12 

125.369 

K2V 

12 

6.992 

K3V/K5V 

12 

1.561 

K7Ve/wdM4 

13 

104.021 

G1III/K0III 

14 

0.814 

dMle/dMle 

15 

3.276 

M2V/M4Ve 

15 

46.817 

K2-3V/K2-3V 

15 

1.268 

M4Ve/M4Ve 

16 

5.975 

K4V/K7.5V 

17 

1.788 

dK5e 

17 

14.300 

K2V/? 

20 

133.286 

G5V 

21 

2.082 

{F7V/JG4V 

21 

. 140 

F6V/G0V 

21 

13.198 

K3-4V-IV 

23 

3.033 

dK2/dK2 

23 

1.034 

dM0e/(dM5) 

23 

20.521 

G8IV-III 

24 

2.879 

K0V/(K5-M2V 

24 

2.170 

dMle/dMle 


31 


198.233 

-59.817 

HR 4980 

25 

4.233 

GOV/GOV 

354.879 

28 . 246 

KT Peg 

25 

6.202 

G5V/K6V 

110.090 

-52.309 

HR 2814 

26 

122.169 

F9.5V{K3:V/(K5V 

286.885 

30 . 254 

V478 Lyr 

26 

2.131 

G8V/(dK-dM) 

74.109 

64.403 

BD+64 487 

27 

44.380 

K2V/K2V 

140.608 

40.201 

BF Lyn 

29 

3.804 

K2V/(dK) 

185 . 549 

73 . 248 

AS Dra 

29 

5.413 

G4V/G9V 

358.767 

28.634 

II Peg 

29 

6.724 

K2-3V-IV 

320.256 

40.345 

HR 8170 

29 

3.243 

F8V/wK5V 

239.683 

25 . 570 

MS Ser 

30 

9.015 

K2V/K6V 

88.363 

-43 . 558 

SZ Pic 

30 

4.960 

G8V 

272.067 

29.691 

V815 Her 

31 

1.810 

G5V/(M1-2V) 

11.835 

24.267 

zeta And 

31 

17.769 

/K1III 

150.007 

24.553 

DH Leo 

32 

1.070 

{K0V/K7V}K5V 

61.884 

-52.568 

AG Dor 

32 

2.562 

KIVp 

159.009 

-11.909 

LR Hya 

34 

6.866 

KOV/KOV 

325.386 

-14.047 

42 Cap 

34 

13.174 

G2IV 

302.112 

15.675 

BD+15 4057 

34 

5.434 

G5V/G5V 

67.155 

19.741 

vB 69 

35 

41 . 660 

KOV 

40.856 

-37.928 

UX For 

35 

0.955 

G5-8V/CG) 

54.197 

0.588 

V711 Tau 

36 

2.838 

G5IV/K1IV 

176.996 

20.219 

93 Leo 

36 

71.690 

A6:V/G5IV-II 

120.649 

57 . 274 

54 Cam 

38 

11.068 

F9IV/G5IV 

259.357 

-66 . 949 

V824 Ara 

39 

1.682 

G5IV/K0V-IV 

273.795 

-32.789 

CP-32 5229 

40 

0.000 

K5Ve 

63.108 

23.675 

BD+23 635 

40 

2.394 

dKO/dMO : 

271.457 

21.447 

ADS 11060C 

42 

25.762 

K7:V/K7V 

271.457 

21.446 

V772 Her 

42 

0.880 

{GOV/ M1V }G5V 

72.298 

24.812 

V808 Tau 

42 

11.929 

K3V/K3V 

65.684 

15.056 

BD+14 690 

45 

4.000 

GOV 

64.412 

16.948 

V818 Tau 

45 

5.609 

G6V/K6V 

64.507 

18.257 

BD+17 703 

46 

75.648 

G4V/G8V 

345.115 

-33.743 

TZ PsA 

46 

1.643 

G5Vp 

315.608 

27.807 

ER Vul 

46 

0.698 

G0V/G5V 

332.170 

45.742 

AR Lac 

47 

1.983 

G2IV/K0IV 

88.085 

-57.156 

HR 2072 

48 

0.000 

F/G5-8III 

51 . 647 

28.715 

UX Ari 

50 

6.438 

G5V/K0IV 

67 . 987 

36.743 

V492 Per 

50 

21 . 290 

K1III 

47.170 

25.592 

BD+25 497 

50 

6.150 

G4V/G6V 

343 . 260 

16.841 

IM Peg 

50 

24.650 

K2III-II 

319.645 

11.569 

BD+10 4514 

50 

3.966 

{F9V/G0VJGIV 

323.050 

0.222 

BD-00 4234 

50 
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